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Abstract 

We obtain C ' a regularity estimates for nonlocal elliptic equations that are not necessarily 
translation invariant using compactness and perturbative methods and our previous regularity 
results for the translation invariant case. 

1 Introduction 

We continue with the program started in [2] to extend the existing theory for (fully non linear) 
second order elliptic equations to nonlocal equations. In this paper we obtain nonlocal Cordes- 
Nircnbcrg type estimates. In [5] we obtained interior C 1,a estimates for nonlocal equations that 
are translation invariant. In this paper we extend those results by perturbative methods to non- 
translation invariant equations and other situations that were not covered in the results in [2] . 

The main theorem in this paper (Theorem 15. 2|) is presented in a somewhat abstract setting 
in order to make it applicable to as many situations as possible. It is a general perturbative 
argument that says that if an equation is uniformly close to another one which has C l a solutions, 
then the solutions of the former are also C 1,a . In the last section of the paper we present concrete 
applications. The equations we have in mind are the purely nonlocal Isaacs equations, which have 
the form 

I(«, x) := inf sup f(u(x + y) + u(x - y) - 2u{x)){2 - er f^^ dy = f(x) in fi. (1.1) 

R" 

where a a p is a family of nonnegative functions with indexes a and (3 ranging in arbitrarily sets A 
and B. This type of integro-differential equations appear in stochastic games with jump processes. 
For different assumptions on the kernels a a p we obtain different regularity results. Some of the 
most interesting applications in the last section of this paper include the following. 

• if < A < a Q/ 3 < A and |Vya a) g| < C\y\ , and a a p(x,y) is continuous in x for a modulus of 
continuity independent of y, then a solution u of (jl.ip is C 1,e for some e > in the interior 
of O (with an estimate depending on ||u|| LOO ). 

• if < A < a a p < A and |Vj,a Q ^| < C| y | — 1 , a Q) g is constant in x, and A — i] < 1 < X + r/ for a 
small enough rj, then the solution u of (jl.l|) is C 2,£ for some e > in the interior of f2 (with 
an estimate depending on |M| LOO ). 

The equations we are most interested in arise from stochastic control problems and stochastic 
games (See [5]). However we develop a general theory that applies as soon as some abstract 
assumptions are satisfied by a nonlocal equation. Note that in the expression (|1.1[) . it only makes 
sense to consider symmetric kernels a a p{x, y) = a a p{x, —y). Indeed, the antisymmetric part of the 
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kernel would have no effect when we integrate it against the second order incremental quotients 
of u (which are inherently symmetric). 

In [T] the regularity of viscosity solutions of nonlinear integro-differential equations that are 
not translation invariant is also studied. They obtain Holder continuity of the solutions. Their 
hypothesis are quite different from ours. They require certain modulus of continuity of the x 
dependence in the kernels. Their equations also involve second and first order terms, and they 
are able to include some degenerate equations into their hypothesis. In [J a parabolic and convex 
integro-differential equation is considered and it is proved that the solutions arc Lipschitz in space 
and C 1 / 2 Holder continuous in time under appropriate conditions. 

In the translation invariant case, the C 1 '" estimate is a somewhat general fact for equations 
satisfying a comparison principle due to the fact that the directional derivatives of solutions satisfy 
some (a priori discontinuous) elliptic equation. If we try to apply the same idea for x-dependent 
elliptic equations we usually need to impose an excessive amount of regularity on the dependence. 
An alternative is to use compactness properties of equations and solutions (see [3] and [I]). In 
the local (classical PDEs) case the approximation takes place in a bounded domain (say the unit 
ball) and simply consists in showing (by compactness) that if the oscillation of the dependence on 
x of the variable equation is small enough, any bounded solution is close to the solution to the 
translation invariant equation inheriting its regularity up to a small error. 

A good portion of the theory depends very little on the actual integral form of the operator 
I(u, x) and it gives rise to interesting applications. Therefore, more generally, by a nonlocal 
operator I we will mean a rule that assigns a value l(u,x) for any function u that is C 2 (x) and 
docs not grow too much at infinity. A safe assumption would be that u is bounded since the kernels 
decay at infinity, but for particular cases it is enough to assume that u £ L 1 (K. n ,uj) for a suitable 
weight uj related to the size of the nonlocal interactions (more concretely, a bound for the tails in 
the kernels K a fj). We also assume that l(u,x) is continuous if u e C 2 . Degenerate ellipticity is 
understood in the following way, if v > u and v{x) = u{x) for some x then I(v,x) > I(u,x). Note 
that the usual elliptic second order operators can be thought of as a special case of this nonlocal 
elliptic operators, and also as a limiting case of integral operators of the form written above. 

For every fixed Xq, the operator 1^,, = I(— ,x<j) given by I Xo (u,.t) = I(t x _ Xo u, xq), where 
T x u(y) = u(x + y), is translation invariant. This is just notation, which corresponds to freezing 
the coefficients in the classical linear theory. 

As mentioned above, we need compactness of the family of solutions to equations with depen- 
dence on x, like De Giorgi or Krylov-Safonov theory for local operators. This we are able to do 
in the case of uniform ellipticity with respect to a family of linear operators. We say that I is 
uniformly elliptic with respect to a family of linear operators C if 

M c v(x) < l(u + v,x)- l(u, x) < M^v(x) (1.2) 

where 

y[~^v(x) := inf Lv(x) 
M~£v(x) := sup Lv(x) 

LeC 

This is the same as the definition in [2J, but now we do not take the operator I necessarily 
translation invariant, and we want to include the case that I may be a second order partial 
differential operator (local operator). For second order equations, this is in fact the definition in 
[5] with respect to the Pucci extremal operators. 

The class C represents the set of kernels that are allowed in an Isaacs equation like (jl.ip . 
Depending on the type of regularity results that we want to obtain, we need to impose less or 
more conditions on these kernels. That is the reason why we need to keep a flexible definition 
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of cllipticity with allows us to consider different classes C A common class of operators is the 
class Cq of uniformly elliptic operators of order a defined in [2] and also later in this paper. These 
are the linear operators whose kernels are comparable to that of the fractional Laplacian of order 
a but no regularity is imposed respect to y. It is the class for which Holder estimates and the 
Harnack inequality is proved in [2J. However, in order to prove C l a estimates for a nonlinear 
problem like a more restrictive class is needed which imposes some regularity of the kernels 

respect to y. 

For the class Cq, or any subclass of it, we can allow the functions u to have certain growth 
at infinity. The exact condition would be that l(u,x) is well defined as long as u G C ' (x) and 
u 6 L 1 (R™, ) meaning that 

/\u(y)\ r~, — i — dty < +00. 

Clearly, assuming that u is bounded is a sufficient condition but it is not necessary. In particular, 
for a > 1, we can still compute l(u, x) if u has linear growth at infinity. 

The main result of this paper is that if some fixed equation l^u(x) = has interior C 1 '" 
estimates, and I(— ,x) stays close to I^ ^ in every scale as x varies in a given domain, then the 
equation l(u, x) = f{x) has interior C 1,a estimates for some a < a. This is a nonlinear fractional- 
order version of the classical Cordcs-Nircnbcrg estimate for second order elliptic equations. 

For linear elliptic PDEs, the classical Cordes-Nirenberg estimate ([5], [5]) says that an equation 

aij(x)diju(x) = f(x) 

has interior C 1,a estimates, provided that / £ L°° and \a,ij(x) — 8ij \ < e for a small enough e. In 
[3], this regularity result was generalized to fully nonlinear equations. Our results in this paper 
recover the second order case since our estimates are uniform with respect to the order of the 
equation. 

The main theorem in this paper is given in section [5J The proof follows more or less the idea of 
the proof in [4] chapter 8, but taking care of the nonlocal effects of the values of the function away 
from the domain of the equation. It is an approximation type argument, based on Lemma 14.71 in 
section |4] that says essentially that if two equations are very close, then so are their solutions. 

It is necessary to give a meaning for two nonlocal equations to be close that will not impose 
restrictive conditions in the long range interactions and will still provide the necessary compact- 
ness. For this purpose, a suitable norm is defined in a preliminary way in section [2] for nonlocal 
operators. In section[5l the definition of the norm is refined so that it represents that two operators 
are close even when we rescale them. 

The proof of Lemma l4.7l is by compactness, so it requires that the equations that we deal with 
have continuous solutions up to the boundary. This is why we have section [3] where we construct 
appropriate barriers to show that solutions to uniformly elliptic nonlocal equations arc continuous 
up to the boundary. 

Our main theorem is presented in a very general form in section [5l So in section [6] we provide 
a list of more concrete applications. The study of the regularity properties of nonlinear integro- 
diffcrcntial equations is a relatively new subject. This paper is part of a program trying to extend 
the fundamental regularity results of fully nonlinear elliptic PDEs to nonlocal equations. The 
family of integro-differential equations is very rich in the sense that we can come up with a very 
large variety of equations. We try to provide theorems that apply to a general class of equations. 
There is always a balance between generality of results and simplicity of exposition and at this 
moment the choice is not always easy to make. 
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2 Preliminaries: Convergence of operators, uniform ellip- 
ticity and viscosity solutions. 

In this section we introduce and discuss closeness of operators, uniform ellipticity and the notion of 
viscosity solutions. We also recall the regularity results of [5] and [5] for uniformly elliptic integro- 
diffcrential equations. We also recall the interior C 1 '" estimates from [5] with the hypothesis 
slightly relaxed. 

We start with a definition of nonlocal operators (Definition 12. We next describe in what 
sense we will say that two operators are close. Since our approach is based on the viscosity method, 
closeness is measured by testing against smooth functions (Definition 12.21) . Next we define and 
discuss in what context our operators are uniformly elliptic. At that point we are ready to define 
what a viscosity solution to our operators is (Definition 12.51) . Then we review the regularity 
theorems proved in [2] and [9] about viscosity solutions. 

In [9] , there is a very simple proof that solutions of integro-differential equations with variable 
order are Holder continuous, as long as the order of the equation stays strictly away from and 
2. The more complicated regularity estimates in [2] hold even when the order of the equation 
gets arbitrarily close to 2, however, only equations with constant order are considered. To a great 
extent, this was a choice for simplicity in the presentation, since many of the results in [5] would 
hold for equations of variable order, but the proofs would be somewhat messier. The results we 
present in this section are only a small variation of the results in [2J. For our main theorem in 
sectional we will make a perturbation argument around an equation with interior C 1,a estimates. 
The applications in section [5] would extend with any extension of our theorems 12.61 and 12 . 71 below. 

We will consider an absolutely continuous weight ui depending on the type of operators that 
we want to deal with. The purpose of this weight is to measure the contribution of the tails of 
the integral in the computation of an integro-differential operator lit. Typically we have in mind 
= 1/(1+ \y\ n+cr ") f° r some u o G (1,2), so that u G L 1 (R", us) means that the integrals in (fTTTI) 
are not singular at infinity for any a > cr . Throughout this paper we always make the following 
assumptions on us. 



These are the only things we need to assume on the weight us for the general theorems of 
section |4] The results in this paper do not apply if the order a of the equation is smaller than 
one, in that respect assumption (|2.1|) is reasonable. Assumption (|2.2|) is more technical and rules 
out singular measures us. 

We use the following convenient definition for a nonlocal operator. 

Definition 2.1. A nonlocal operator I is a rule that assigns to a function u a value I(u,x) at 
every point x, satisfying the following assumptions. 

• I(u,x) is well defined as long as u G C 2 (x) and u G L 1 (R",ij). 

• If u G C' 2 (il) n L (M. n ,ui), then I(u,x) is continuous in il as a function of x. 

By u G C 2 (x) we mean that there is a quadratic polynomial q such that u(y) = q(y)+o(\x — y\ 2 ) 
for y close to x. 

Definition 12.11 is very general. It becomes useful when complemented with the concept of 
(nonlocal) uniform ellipticity (jl.2[) . We do not include (|1.2I) in the definition above because we 
want to define the norm of a general operator (defined below) and compute it for differences of 
elliptic operators. 



i + Hgl^R",^) 

sup w < C{r)w{y) 
B r (y) 



(2.1) 
(2.2) 
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The results in this paper are based on approximation techniques. In order to use that type of 
argument, we need to understand when two nonlocal operators are close to each other. Since we 
use the Crandall-Lions viscosity concept, it is natural to define closeness of operators by testing 
them against the appropriate family of functions. We define the following norm. 

Definition 2.2. Given a nonlocal operator I, we define ||I|| in some domain £1 with respect to 
some weight lu as 

||I|| := sup {|I(u,ar) |/(1 + M) : ieS] 

u e C 2 (x) 
IMIlmr-^) < M 

\u(y) — u(x) — (y — x) ■ Vw(x)| < M\x — y\ 2 for any y £ B\(x)} . 

The idea of this definition of ||I|| is that j|Ifc — I| — > if Ik(u, x) — > I(u,x) in a somewhat 
uniform way. The condition ||u|| L i^ R „ w % < M could be replaced by |u| < M in R n which would 
make convergence weaker. 

Note that in the local case when I(u, x) = F(D 2 u, x), we have w = and ||I|| = (3(F), where (3 
is the one from 0]. 

Our uniform ellipticity assumption (|1.2[l depends on the class £ of intcgro-diffcrcntial operators. 
A class £ is a set of operators L of the form 

Lu(x) — J (u(x + y) + u(x — y) — 2u(x))K(y) dy, 

for a family of nonnegative kernels K satisfying J min(l, \y\ 2 )K(y) dy < +oo. 

We will make some minimal assumptions on every class C. In this whole paper, every class C 
will satisfy the following two assumptions, that connect the class with the choice of the weight uj. 

Assumption 2.3. If K = sup K a is the supremum of all kernels in the class, then for every r > 

K(y) < CMv) *f \V\ > r. 
for some constant C r (depending on r). 

Assumption 2.4. For some constant C and every L G C, we have \\L\\ < C (which clearly implies 
||M+|| < C and also ||M~|| < C). 

Assumption 12 . 31 savs that the weight u controls the tails of the kernels at every scale. Assump- 
tion 12.41 says that near the origin the kernels are uniformly controlled by C 2 functions. These 
are very weak assumptions and are all that is needed for our stability results of section UJ Our 
regularity results, however, require stronger assumptions in the class £. 

In order to understand the assumptions better, let us consider C to be the set of all fractional 
Laplacians — (— A) CT / 2 for a £ [ao, 2], it would satisfy both assumptions with U)(y) = (2 — co)/(l + 
|y| n+<7 °). On the other hand if we take the set C of all operators with kernels of the form l/\y\ n+a 
(where we arc omitting the normalizing factor « (2 — cr)), then this class would not satisfy 
Assumption 12.41 

In the case of local equations, the class £ would consist of linear second order partial differential 
operators, and the weight u) is irrelevant (it can be taken as cj = and also K = 0). Our results 
in this case reduce to the usual results for uniformly elliptic partial differential equations [3]. 

We recall the definition of viscosity solutions given in [2]. 
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Definition 2.5 (Viscosity Solutions). Let I be a uniformly elliptic operator in the sense of (|1.2[) 
with respect to some class C. A function u : R ra — >• M., upper (lower) semi continuous in O, is said 
to be a subsolution (supersolution) to I(u,x) = f{x), and we write l(u, x) > f(x) (l(u,x) < f{x)) 
in fl, if every time all the following happens 

• xq is any point in D,. 

• N is a neighborhood of xq in Q. 

• ip is some C 2 function in N. 

• (p(x ) = u(xo). 

• f(x) > u(x) (<f{x) < u{x) ) for every x G N \ {xq}. 
Then if we choose as test function v, 

{ip in N 
u m I"\JV, 

we have l(v,x ) > f(x ) (l(v,x ) < f(x ))- 

A solution is a function u which is both a subsolution and a supersolution. 

The first theorem we borrow from [2] is the interior Holder estimates. We change the assump- 
tion of u being bounded for u € L 1 (M. n , 1/(1 + \y\ n+aa ))- Let us first recall the class Co of integro 
differential operators. 

An integro differential operator L belongs to Co if its kernel K satisfies 

(2-a)— ^— < K(y) < {2-a)—^— . (2.3) 
In this case and take a very simple form: 
M+ v(x) = (2 

M C( v(x) = (2 

R" 

where 5 represents the second order incremental quotients 5u(x, y) = u(x + y) + u(x — y) — 2u(x). 

The class Co depends on the parameter a G (0,2). Whenever it is necessary, we will stress 
that dependence by writing Co (a). 

The next theorem concerns equations with kernels that can be arbitrarily discontinuous in 
terms of x. In particular, it applies to the difference of two solutions of a nonlinear integral 
operator. 

Theorem 2.6. Let a > a for some a > 0. Let u be function in L 1 (R n ,u) forcj = l/(l + \y\ n+a "), 
such that u is continuous in B\ and 

M^ o u > -C in Bi 
M £n u < Co in Bt 

then there is an a > (depending only on X, A, n and o~o) such that u G C a (Bi/ 2 ) and 

«C«(B 1/2 ) < C( sup \u\ + |M| L i (R „ jtlj) + Co) 
Bi 

for some constant C > which depends on o~o, A, A and dimension, but not on a. 



ASv(x,y) + — X5v(x,y)~ 
\y\ n+a 

X8v(x, y) + — A6v(x, y)~ 



dy 



dy 



(2.4) 
(2.5) 
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Proof. Consider u(x) — xb 1 {x)u{x). A direct computation shows that M + w > — Co — c |M| L i(- R „ ^ 
and M~m < Co + c |M| L i( K „ w ^ in -B3/4. We can now apply theorem 12.1 from [2] to u. □ 



Thcorcm l2.6l is a nonlocal version of Krylov-Safonov theorem for elliptic PDEs in nondivcrgence 
form. A nonlocal version of De Giorgi-Nash-Moser theorem can be found in [5] . 

We will also mention, as a consequence of the previous theorem, the interior C 1,a estimates 
from [5] for functions with some growth at infinity. 

Given po > 0, we define the class £* by the operators L with kernels K such that 

\VK{y)\<Cu(y) inM"\i? P0 (2.7) 

Theorem 2.7. Assume a > ctq > 0. There is a po > (depending on X, A, (To and n) so that 
if I is a nonlocal translation invariant uniformly elliptic operator with respect to C* and u is a 
continuous function in B\ such that u £ L 1 (IR™,a;) for cj = 1/(1 + lu = in B\, then 

there is a universal (depends only on X, A, n and o~$) a > such that u £ C 1+a (B 1 / 2 ) and 



u c^(B 1/2 ) < C (jmp |u| + \\u\\ L i [mnM) + \I0\j 



for some constant C > (where by 10 we mean the value we obtain when we apply I to the function 
that is constant equal to zero). The constant C depends on X, A, a®, n and the constant in (|2.7[) . 



Sketch of the proof. The proof of this theorem goes along the lines of the proof of Theorem 12.1 
in [2] by applying Theorem 12.61 to the differential quotients w h = (u(x + h) — u{x))/\h\P for 
j3 = a, 2a, . . . , 1. In each step, the quotient w h is written as w h = w'^ + w\ where w\ vanishes in 
a ball B r (with 1/2 < r < 1) and w*l is bounded in R™, M^ 1 < C and M+w'l > -C in B r . This 
constant C, like in the proof of theorem 12.1 in [2], comes from the estimate 

. . ,.\K(y) - K(y - h)\ , 

Hx + y)\ i ii dy < C. 

\h\ 

In [2] the assumption was u £ L°° whereas now we have u £ L 1 (R n ,w) and so is why the 
assumption (12.2) in [2] necessary for the bound of the constant C is different from our assumption 
(|2.7[) in this paper. □ 

Note that sup Bi \u\ + \\u\\ L1 ^ Rn ^ < Csup |n(ar)|/(l + |x| 1+Q ) for any a < (Jo- 
in comparison with Cq, the class £* adds the extra regularity assumption (|2.7[) to the kernels 
with respect to y. This condition is needed in order to prove the interior C 1,a estimates in the 
ball of radius one. However, the hypothesis is not scale invariant, which makes it not suitable to 
the type perturbation arguments used later in this paper. Therefore, in section [5] we will further 
restrict this class in order to make it scale invariant. 

Remark 2.8. Assumption (|2.7[) is used to control the effect on the regularity of u at x of the 
points y that are at least po away from x. This condition depends on the type of control we have 
of u far from x. In [2] we were assuming the stronger condition u £ L°°(M. n ), so we had the weaker 
condition for the kernels 



\K(y)-K(y-h)\ p 

j— I dy < C every time \n\ < — 

\n\ 2 



PO 
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Now we have the weaker condition u G L 1 (R",o;) instead. If we want our estimate to depend on 
different norms of u, our condition (|2.7|) will change accordingly. 

It is also convenient to define and as 

Mtu(x) = lim Mf , Mx) 

M^u(x) = lim M7 , \U(x) 
* v ' <r->2 £ °(°") 

A simple computation shows that Af^~u(x) is a second order uniformly elliptic operator of the 
form M^u(x) = F{D 2 u) with 

F(A):= J A(y,Ay) + -X(y,Ay)- AS. 

Note that the ellipticity constants A and A of as a partial differential operator depend only on 
A, A and dimension. Also F(0) = 0, so A/ 2 + u(x) < M + (D 2 u) where M + (D 2 u) is the usual Pucci 
operator with ellipticity constants A and A. Naturally, we also have the corresponding relations 
for M 2 ". 

3 Boundary regularity 

In this section we show that if we have a modulus of continuity on the boundary of the domain of 
some equation, then we can find a (possibly different) modulus of continuity inside the domain. 
This is done, as usual, by controlling the growth of u away from its boundary data through barriers, 
scaling and interior regularity. 

We start by showing that a truncation of an appropriate radial function is a barrier outside 
the unit ball. Our barrier function is good as a supersolution M + <^ < for all values of a greater 
than a given Co, where M + is the maximal operator Mj^^. Another way to say this would be to 

define a larger class C which is the union of all classes Co(a) for a G (ctq, 2), then M.~}.(p < 0. 

Using this barrier function we will prove that the solutions of uniformly elliptic equations are 
continuous on the boundary, and then using their interior Holder regularity (Theorem I2.6|) we 
get that they are continuous up to the boundary in the domain of the equation. The result in 
Theorem l3.3l is restricted to equations of constant order only because Theorem l2.6l is. In fact, the 
same continuity result holds for operators of variable order as long as the order remains strictly 
in between and 2 (because of the results in [3]). 

Lemma 3.1. Given any <7q G (0,2), there is an a > and r > small so that the function 
u{x) = ((\x\ — 1) + ) Q satisfies AI + u(x) < in B\ +r \ B\ for any a > o~q. 

We defer the proof of this lemma to the appendix. We use it to obtain the following corollary. 

Corollary 3.2. For any constant C, there is a continuous function ip such that 

• tp = in B\ . 

• tp > in R". 

• cp > 1 mR n \ B 2 . 

• M+ip < in R" \ B 1 for any a > a - 
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Proof. Consider ip = min(l, C((|a;| — l) + ) a ) for some large constant C and apply Lemma |3~T1 □ 

We will use the function ip of Corollary 13.21 as a barrier to prove the boundary continuity of 
solutions to nonlocal equations. We will do it in spherical domains. Note that cp is a supersolution 
outside of a ball, so the same method could be used in domains with the exterior ball condition. 

Theorem 3.3. Let a > tro > 0. Let p be a modulus of continuity and M + and M~ be the maximal 
operators M^T and M^^. Let u be a bounded function such that 

M+u{x) >-C inB x 
M~u(x) <C inBx 
\u(y) — u[x)\ < p(\x — y\) for every x € dB\ and y G R™ \ B±. 

Then there is another modulus of continuity p so that \u(y) — u(x)\ < p(\x — y\) for every x G B\ 
andy G R". 

The modulus of continuity p depends only on p, X, A, o~o, n, \\u\\ LX> and C (but not on a). 

The modulus of continuity p in Theorem 13.31 is in some sense the worst between the internal 
regularity and the boundary regularity. 

Before proving Theorem 13.31 we mention a corollary. 

Corollary 3.4. Let a > do > 0. Let p be a modulus of continuity. Let u be a function such that 

IMlL 1 (i?",l/(l+|'y|" + "o)) - C 

M + u{x) > -C in B Y 

M~u(x) <C mBx 

\u(y) — u{x)\ < p{\x — y\) for every x G dB\ and y G K™ \ B\. 

Then there is another modulus of continuity p so that \u(y) ~ u(x)\ < p(\x — y\) for every x G B\ 
and y G R™. 

The modulus of continuity p depends only on p, A, A, o~q, n, and C . 

Proof of corollary. Since \u{y) — u(x)\ < p(\x — y\) if x G dB\ and J B \u\ dx < C, then |m| < C 
on dB\. Thus, from the modulus of continuity, u < C in i?2- 

We apply Theorem [3731 to u(x) = min(C, max(u, C)). Note that M+ii < M+u + M+(u - u) 
and M + (u - u) is controlled by Hull^^n ,i/(i+|j,|«+<'o))- D 

In order to prove Theorem 13. 3[ we will need a couple of preparatory lemmas that mirror the 
second order approach. Lemma 13.51 reconstructs the modulus of continuity on the boundary, and 
Lemma l3.6l in the interior. 

Lemma 3.5. Assume a > o~q > 0. Let p be a modulus of continuity. Let u be a bounded function 
such that 

M+u(x) > -C 

u(y) — u{x) < p(\x — y\) for every x G dBi and y G R™ \ B\. 

Then there is another modulus of continuity p so that u{y) — u{x) < pQx — y\) for every x G dB± 
andy G R". 

The modulus of continuity p depends only on p, X, K, ao, n, the constant C above, and \\u\\ Lao . 
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Proof. We first observe that since a > ao the function p(x) = max(0, 4 — |x| 2 ) satisfies M + p > —c 
in B\ for some c > 0. So we can consider u — ^p and reduce the problem to the case C = (the 
values of |M| ioo and p are affected depending on C). 

Let po be the modulus of continuity of the function ip of Corollary 13.21 Let p be 



p(r) = inf (p(3fl) + ||u|| LOO Po (£) 



Let us first show that p is indeed a modulus of continuity. The function p is clearly monotone 
increasing because po is. We still need to show that for every e > 0, there is an r > such that 
p(r) < e. We first choose R such that p(3R) < e/2, and then choose r such that ||u|| i0 o Po (jj) < 
e/2. 

Now, we show that p(\y — x|) > u(y) — u(x) for any y £ M™ and x e For any i? > 0, 

Consider the barrier function 

B(y) = u(x) + p(3R) + \\u\\ Lae ^ [-x + V -^- 

We see that B(x) > u{x)+p(3R) > u in B 3R (x)n{W n \B 1 ). Also B(x) > \\u\\ Laa in R n \B 3R (x). 
Moreover M + B < 0, whereas M + u > in B\, so u < B < u(a;) + — x\) everywhere. Taking 
infimum from all possible choices of R we obtain p(\y — x\) > u(y) — u(x). □ 

Lemma 3.6. Assume a > o~q > 0. Let p be a modulus of continuity. Let u be a bounded function 
such that 

M + u(x) > -C in B 1 (3.1) 

M~u{x) <C inB x (3.2) 

\u(y) — u(x)\ < p(\x — y\) for every x G dB\ and y 6 R™. (3-3) 

XTien i/iere is another modulus of continuity p so that \u(y) — u(x)\ < p(\x — y\) for every x £ B\ 
and y£R n . 

The modulus of continuity p depends only on p, X, A, the constant C above, ao, n, and \\u\\ L<x> . 

Proof. Let x £ B\ and y G W" . We will consider two cases. If dist(x, y) < dist(x, dB\)/4 or if 
dist(x,y) > dist(x,a_Bi)/4. 

If dist(x, y) > dist(x, dBi)/A then let xo be a point in dB\ where dist(x, Xo) = dist(x, dBi), 

\u(x) - u(y)\ < \u(x) - u(x ) \ + \u(x ) - u(y)\ 
< p(\x - x \) + p(\x - y\) 
<p{\x-y\/A) + p{S\x-y\/A) 

Let us consider now the case dist(a;, y) < dist(x, 9J5i)/4. Let r = dist(x, dBi)/2. Let us 
truncate the function u at u(xo) ±p(4r) by defining u := min(M(xo) +p(4r), max(u, u(xq) — p(4r))). 
So we have 

\u{z) -u(z)\ < min((/o(r + |z - x|) - p(4r)) + ,M) 

where M = 2 \\u\\ Lao . 

We now scale the functions so that B r (x) becomes B\. We write 

v(z) = u(x + rz) 
v(z) = u(x + rz) 
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and we have 

\v(z) -v(z)\ < min((p(r + r\z\) - p(4r))+,M). 

We want to estimate M + v in B\. We know that M + v > —Cr a in B\ by rescaling (|3.1[) . Also 
M + v > M + v — M + (v — v). We now estimate M + (v — v) in B\ by above. We have 

tl/j+1 t \ -f^^ f0 \ a /" 2min((p(r + r|y|)-p(4r))+ M) 
A/ + (u(z) < (2-cr)A J dy 

=:I r 

We observe that I r — > as r — > uniformly in cr as long as cr > cr . We thus have —Cr a —I r < M + v 
in Bi, and Cr a + I r -> as r -> 0. 

We can repeat the same computation to obtain < Cr' 7 + I r in Si. Applying interior C a 

estimates [2], we obtain 

v(z) - TJ(0) < C(Cr a +I r + p(4r))\z\ a = Cm r \z\ a 

where m r := Cr a + I r + p(4r) goes to zero as r — > uniformly for cr > oq. 

Now we scale back with y — x + rz to apply the above estimate to u. We have 



\u{y) — u{x)\ < Cm r 



\x - y\ 



so we have proven that \u(y) — u{x)\ < p(\x — y\) as long as wc can chose p(d) > C swp r>id d ™ r . 
There will be such modulus of continuity as long as we can prove that 

sup -> as d -> (3.4) 

r>4d ^ Q 

We will prove this by contradiction to finish the proof of this lemma. Assume (|3.4[) was not 
true. Then there would be an e > and two sequences Vi and di with r*j > 4di and di — > such 
that dfm r Jrf > e. Since m r is bounded above, df/rf cannot go to zero. But that would imply 
that Vi —> and dfm r Jrf < m r jA a — > 0, which gives a contradiction. □ 

Remark 3.7. If we only considered Holder modulus of continuity, the proof would not be much 
shorter but it would be more direct. It is nevertheless interesting that it can be made to work for 
arbitrary modulus of continuity. 



Proof of Theorem \3.S\ We apply Lemma l3.5l to both u and — u to obtain a modulus of continuity 
that applies from any point on the boundary dB\ to any point in space R™. Then we use Lemma 
13.61 to finish the proof. □ 

Remark 3.8. The conditions M + u > —C and M~u < C are implied by u being a solution of a 
linear equation like 

J 6u(x, y)(2 - <7)^p^ dy = f(x) 

R" 

no continuity whatsoever is required for the coefficient a(x,y), just that A < a(x,y) < A. The 
function / is only required to be bounded. The modulus of continuity is independent of a as long 
as a > o"o > 0, but the result was proved for constant order a. We could obtain a corresponding 
result for a solution u of an equation with variable order o~{x) if instead of Theorem 12 .61 we used 
the estimates from [9]. The disadvantage is that those estimates would blow up as a —> 2. 
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4 Approximation results. 



The purpose of this section is to show that if two equations are very close to each other, then so 
are their solutions. The proof is by compactness: if a sequence of solutions and corresponding 
equations converges, the limiting function is the unique viscosity solution of the limiting equation. 
We stress that, thanks to the concept of viscosity solutions, it suffices that the convergence of the 
functionals be weak (i.e. against smooth test functions). 

We obtain our desired result in great generality in Lemmas 14. 31 and !4.7l To prove those lemmas 
we need some previous technical results that we develop in this section. 

In the following lemma, we notice that the test functions used in the definition of viscosity 
solutions can be reduced to only the polynomials of degree two. 

Lemma 4.1. In Definition \2.5\ it is enough to consider test functions ip which are quadratic 
polynomials and neighborhoods N which are balls centered at the touching point. 

Proof. Assume Definition ^. 5l holds every time (f is a quadratic polynomial and N is a ball centered 
at the touching point. We will show it also holds for any C 2 function ip. 

Let ip be any C 2 test function from above in a neighborhood N (similar reasoning would apply 
for test functions from below). Without loss of generality we assume ip touches u at the origin. 
Let p(x) be the following quadratic polynomial. 



where in the last inequality we used that ||Mj|| < C. 

Se we have 1(^,0) > l(v p , 0) - Ce > /(0) - Ce for arbitrary e > 0. Thus l(v v ,0) > /(0). □ 

We want to show that if Ifc(itfc, x) = fk(%) and I& — > I, Ufc — > u and //. — > f in some appropriate 
way, then l(u, x) = f{x). It is indeed true that the convergence I& — > I can be taken with respect 
to the norm defined in section [2j however we need a stronger result with respect to a weaker 
version of convergence that we define below. 

Definition 4.2. We say that 1^ — > I weakly in if for every xq € f2 and for every function v of 
the form 



where p is a polynomial of degree two and u £ L 1 (M' l ,o;) ; we have lk(v,x) —> l(v,x) uniformly in 
B p / 2 (xo). 



p{x) = ((D 2 (p(0) + el)x, x) + x ■ Zty(0) + <p(0) 



So, p > tp > u in a neighborhood of the origin B r . Let 




We choose r small so that v p (x) < v v (x) + 2e\x\ 2 XB r - 
By assumption, l(v p ,0) > /(0). On the other hand 



I(iV,0) -l(v p ,0) > M c (v ip ~v p ,0) 

> -M+{{v v -v p )-,0) 

> -CM+(2e\x\ 2 XB r ,Q) > -Ce 
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The following lemma is an improvement of Lemma 4.5 from [2J. 

Lemma 4.3. Let \ k be a sequence of uniformly elliptic operators with respect to some class C. 
Let us assume that Assumvtion \2.3\ holds. Let u k be a sequence of lower semicontinuous functions 
in fl such that 

• l k (u k ,x) < f k {x) in tt. 

• Uk — > u in the V sense in SI. 

• Uk —> u in L 1 (M. n ,u>). 

• fk —> f locally uniformly in Q. 

• Ifc — > I weakly in fl ( with respect to uj) 

• \u k (x)\ < C for every x S fi. 
Then lu < f in Q. 

Proof. Let p be a polynomial of degree 2 touching u from below at a point a; in a neighborhood 

B r (x). 

Since Uk L-converges to u in fi, for large k, we can find Xk and dk such that p + dk touches Uk 
at Xk- Moreover Xk — > x and dk — > as k — >• +oo. 
Since l k (u k ,x) < f k (x), if we let 



Vk 



p + d k in B r (x) 

u k mW l \B r {x) 



we have lk{v k ,x k ) < f k (x k ). Clearly v = lim fe ^ +00 v k . 
Let z 6 B r / A (x). We have 

\l k {v k ,z)-l{v,z)\ < \lk(v k ,z)-lk{v,z)\ + \I k (v,z) -l(v,z)\ 

< max (\M+(v k - v)(z)\, \M+(v - v k )(z)\) + \l k (v, z) - l(v, z)\ 

< sup \L(vk - v)(z)\ + |Ifc(v, z) - l{v, z)\ 

L<£C 



( 



< 



\S(v k -v){z,y)\K(y) dy \ + \h(v,z) - I(v,z)\ 

\B r/ 2 




<C\ / 2\v k (y) -v(y)\ sup w(y + z) dy | + C\v k (z) - v(z)\ + \I k (v, z) - I(v, z)\ 

J z£B r/ i 
l n \B r/2 

< C\\v k - w|| L i (R „^) +C\v k (z) -v(z)\ + \h(v,z) - l(v,z) 



Where we used Assumption 12.31 and 

The term |Ifc(u, z) — I(v, z)\ goes to zero uniformly for z G B r / A (x) since we assumed I k — > I 
weakly. The sequence \\v k — v\\ L1 ^ — > and u k — > u uniformly in B r (x), therefore the right hand 
side in the above inequality goes to zero uniformly for z £ -B r /4- We obtain l k (v k ,x) — > x) 
uniformly in B r / A {x). 

We have that Iv is continuous in B r (x). We now compute 

\h(vk,x k ) - I(v,x)\ < \l k (v k ,x k ) - I(v,x k )\ + \l(v,x k ) -l(v,x)\ . 

So l k (v k , x k ) converges to x), as k — > +oo. Since Xk — > x and f k —tf locally uniformly, we 
also have f k (x k ) — > f(x), which finally implies l(v,x) < f(x). □ 
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Remark 4.4. Note that the weight oj is important. Consider the following example. Let 1^ = 
— (— A) 1_1 / fe and I = A. Then I is a local operator, and we may think that we can take u> = 
and overlook the L 1 (R™, cu) convergence. However we can consider the sequence of functions 

u k (x) = Xb 1 (x)(\x\ 2 - 1) - M k XM"\B 2k ■ 

If we choose Mk large enough, then we can make I^Ufe < in B\ for every k, but 1/ = A/ = 2n > 
in Si. 

The point is that in this example there is no single weight uj for which at the same time Uk — > u 
in L 1 (R",a;) and ||I„ — I|| — > with respect to to. 



Lemma 4.5. Let v be a function 

v{x) 



p(x) if x e B r 
u{x) if x ^ B r 



and let L. be a sequence of uniformly elliptic operators respect to some class C satisfying As- 
sumption \2.S[ There is a subsequence Ife. such that fk j (x) := 1^. (v, x) converges uniformly in 

B r /2- 

Proof. All wc need to do is find a uniform modulus of continuity for fk in B r /2 so that the Lemma 
follows by Arzela-Ascoli. 

Recall the notation t z u(x) = u(x + z). Given x, y G B r /2 with \x — y\ < r/8, we have 

fk{x) - fk(y) = h(v,x) - h(v,y) 

< Mj(f - Ty- X V, X) 

Note that v — Ty— X v is just linear function in B r /i(x) 
S(v — T y - X v, x, z)K(z) dz 



< 

R"\B r/4 

where K is the kernel of Assumption 12. 31 thus we have 



< C / (v(x + z) + v(x — z) — 2v(x) — v(y + z) — v(y — z) + 2v(y))uj(z) dz 

K"\B r/4 

< C(\v(x) -v(y)\ + \\T x v~T y v\\ LliRn ^ } ) 



< C J sup \v(x) - v(y)\ + / \t( x - V )V - v\( sup w) dz ] < c(|i - y\ 

\%~~y\<\ x — y\i .r.iyF /i r /2 B r (z) 

Where we use (|2.2|) and c(p) is defined as 

c (p) : = C SU P K^) ~ v (y)l + / \ T (x- y ) v ~ v\u(z) dz 

Clearly, c(p) is a modulus of continuity that depends on v but not on 1^. So the functions 
fk have a uniform modulus of continuity and there is a subsequence that converges uniformly by 
Arzela-Ascoli. □ 
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Theorem 4.6. Let I& be a sequence of operators uniformly elliptic with respect to some class C 
satisfying Assumptions\2.4\ and\2.3[ There is a subsequence It. that converges weakly. 



Proof. In order to make this proof, first we will construct a dense subset of the test functions v of 
the form 

J p(x) if x G B r 
I m(xj it x ^£ is,.. 

This dense subset is constructed as follows. We know that L 1 (M. n ,Lu) is separable, so it has a 
countable dense subset {ui}. Quadratic polynomials are a finite dimensional space which also has 
a countable dense subset {pi}. Now, for every positive integer fc, we construct the functions 



Vk,ii j*2 (*^) 



if x G 5 2 - fc 
zti 2 (a;) if x ^ -B2- & • 



So, given any function v as in (|4.ip and any e > we first choose k such that 2~ fc < r < 2 _fc+1 , 
then we choose i\ and 12 such that ||uj 2 — w|| L i( R „ ^ < £, l-D 2 ^^ — -D 2 p| < £, — < £ in 

B 2 -k, and jp^ — p\ < e in B 2 -k. 

Since the set {i>k,h,i 2 } is countable, we can arrange it in a sequence Vi of the form 



Vi(x) 



Pi (x) if x G B u 
Ui(x) if x (fc B Ti 



such that for every v as in (|4.1[) . there is a such that 

ll«-«i||£i(H»^,) < e 

<r F in R /„ 

(4.2) 



\v — Vi\ < e in B r / 2 



Dv — Dvi\ < e in B r i 2 

r/2- 



\D 2 v-D 2 vA<e ixiB 



By Lemma l4.5[ for each i>j, we can take a subsequence Ifc. such that 1^ converges uniformly 
in B rj / 2 . 

By a classical diagonal argument, there is a subsequence K. such that for every u,, Ife, 
converges uniformly in B r ./ 2 . We call this limit I^t^x). 

If u is any test function, there is an i such that v is very close to Vi in the sense of (|4.2I) . For 
any x G B r / 2 , we have 

Ife(w,x) - l k (vi,x) < M+(w - Wj,x) 

< M+((v - Vi) X B r/2 ,x) + M+((v - Vi )(l - XB r/2 ),x) 

Note that v — Vi is a quadratic function in B r / 2 (x) 

<e||Ife||+ J 8(v - Vi,x,z)K(z) dz 

R n \B r/2 



where K is the kernel of Assumption 12.31 thus we have 

<£||I fe ||+C||w-u ? .|| il(R „ j 
< Ce 
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So, by choosing Vi, we can make Ifc(i>, x) — lk(vi,x) as small as we wish in B r / 2 , uniformly in 
k. So, for j large, 11^ (v, x) — Ioo(«i, x)\ < 2e and thus (v, x) is a Cauchy sequence in L°°{B r / 2 ). 
We define l^v^x) to be the uniform limit of this sequence in B r / 2 . 

Therefore we have shown that Ik{v,x) converges uniformly to IooC^j x) in B r / 2 . To finish the 
proof, all we must show is that this operator loo can be extended to a uniformly elliptic operator 
for all test functions tp. The key is to note that for any two test functions v\ and v 2 of the form 
(|4.1[) we have M^(wi — v 2 ,x) < — v 2l x) < Mj(wi — v 2 ,x), and so this inequality passes to 

the limit and becomes M^(t>i — v 2 ,x) < Ioo(^i — v 2 ,x) < M~^(vi — v 2 ,x). Now, approximating an 
arbitrary test function tp as in the proof of Lemma |4. II we get that there is a unique way to extend 
loo to all test functions <p such that loo is uniformly elliptic respect to C. □ 

Lemma 4.7. For some a > o~o > 1 and a < <jq — I we consider nonlocal operators Io, Ii and I 2 
uniformly elliptic respect to Co(a). Assume also that the boundary problem 

IqU = in Bi 

u = g in E" \ B 1 

does not have more than one solution u for any g continuous and \g(x)\ < M (\x\ + 1) 1+Q . 

Given M > 0, a modulus of continuity p and e > 0, there is an n > (small) and a R > 
(large) so that if u, v, Io, T and I 2 satisfy 

Io(w, x) = in B\, 

Ii(u, x) > —r] in B±, 

l 2 (u, x) < i] in B\ 1 

u = v inR n \B u 

Pi - Io|| < i] in Bi, 

||I 2 -I ||<77 inB u 

\u(x) — u(y)\ < p(\x — y\) for every x £ Br \ B\ and y £ W 1 \ B\, 
\u{x)\ <M(|x| + l) 1+Q . 

then \u — v\ < e in B\. 

Proof. Assume the result was false. There there would be a sequence Rk, Iq , I x , I^', rjk, Uk, 
Vki Ife) fk such that Rk —> oo, rjk — > 0, all the assumptions of the lemma are valid but nonetheless 
sup \u k - v k \ > e in B x . 

Since I fc ' is a sequence of uniformly elliptic operators, by Lemma [4.6l we can take a subsequence 
that converges weakly to some nonlocal operator Iq that is also elliptic respect to the same class 



T (*0 t (*0 



and 



T (fe) T (fc) 



0, then Ij*' and 1^ also converge weakly 



C. Moreover, since 

tO Iq. 

Since both Uk and Vk have a modulus of continuity p on dB\ , then by Lemma 13.61 there is 
a modulus of continuity p that extends to the full unit ball B±. So Uk and Vk have a uniform 
modulus of continuity on Bn k with Rk — > oo. We can find a subsequence that converges uniformly 
on compact sets, thus converging almost everywhere in R n . By dominated convergence theorem, 
Uk and Vk converge in L 1 (R™,w) to some function g. 

Let u and v be the limits of Uk and Vk respectively. Since Uk and Vk converge to u and v 
uniformly in B\ and sup \uk — i>fc| > £ in Si, then u and v must be different. We use Lemma |4. 3 
and we get that u and v solve the same equation Io(m, x) = Io(w, x) = 0. But then u = v which is 
a contradiction. □ 
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Remark 4.8. We will apply Lemma 14.71 to Io translation invariant in which case the uniqueness 
assumption for the boundary problem was proved in [2]. 



We also present the following simplified version of Lemma 14.71 which will not be used in this 
paper but it is interesting by itself. The difference with respect to Lemma 14.71 is that in the 
following version we fix the boundary value g and in exchange we do not require a modulus of 
continuity in Br\ B\ but only on dB\. 

Lemma 4.9. For some a > Co > 1 and a < <To — 1 we consider nonlocal operators Iq, Ii and I2 
uniformly elliptic respect to Cq{o~). Assume also that the boundary problem 

Iqu = in B\ 
u = g inR n \ B 1 

does not have more than one solution u for some given boundary condition g. 

Assume g is continuous on dB\. Given M > and £ > 0, there is an rj > (small) so that if 
u, v, Iq, Ii and I2 satisfy 

Iq{v, x) = in Bi, 

x) > —rj in B\, 

l2(it, x) < 7] in B\, 

u = v = g inl"\Bi, 

||Ii-I ||<7? inBx, 

||I 2 -Io||<»7 inBx, 

then \u — v\ < e in B\ (in this case rj depends on g). 

Proof. The proof follows the same line as the one of Lemma 14.71 We do it by contradiction 
assuming the result was false. There there would be a sequence Rk, Iq , 1^, I2 , ?/fc, Uk, Vk, Ifc, 
fk such that r/k — > 0, all the assumptions of the lemma are valid but nonetheless sup \uk — Vk\ > s 
in Bx- 

The functions and Vk have a fixed value g outside Bx- Since g is continuous on dB\, 
by Lemma 13.61 Uk and Vk are equicontinuous in Bx ■ So there is a subsequence that converges 
uniformly in B\. 

We continue as in the proof of Lemma 14.71 We can take a subsequence such that I^ and 
Iq converge weakly to Io. Let u and v the limits of Uk and Vk respectively, so we know that 
sup \uk — Vk\ > £ ha Bx- But by Lemma 14.31 u and v must solve the same equation Io(it,2;) = 
Iq(v,x) = 0. But then u = v which is a contradiction. □ 



5 C a regularity for variable coefficient equations 

The main purpose of this section is to prove a C 1 '" estimate for nonlocal equations that are not 
necessarily translation invariant. We will give the result in its most general form, and in the next 
section we will provide a variety of applications. 

It is important to stress the concept of scaling. Since our proofs will involve successive rescaling, 
we will need some kind of scale invariance. It is not necessary to require a particular equation to 
be scale invariant, but instead we will consider our equations within a whole class of equations 
that is scale invariant for which our regularity result up to the boundary is supposed to apply, 
and we will jump freely from equation to equation in our rescaling process. We say that the class 
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C has scale a if every time the integro-differcntial operator with kernel K(y) is in C, then the one 
with kernel X n+a 'K(Xy) is also in L for any A < 1. For example the class Cq defined in (|2 . 3[) has 
scale a, whereas the class £* defined in (|2.6p ~ (|2.7p does not because of condition (|2.7p . 

It is a straight forward computation to check that if C has scale a and u solves an equation 
I(u, x) = f(x) that is elliptic with respect to C, then the function w(x) = fiu{Xx) solves a uniformly 
elliptic equation with respect to the same class C, l^^^WyX) = A CT (if(Xx). For example, if I(u,x) 
corresponds to an integro-differcntial operator 



I(it, x) = J (u(x + y) + u(x — y) — 2u(x))K(y) dy 

then I^x corresponds to the integro-diffcrential operator 

I MiA (tt, x) = J{u{x + y) + u{x -y)- 2u{x))X n+ ° K{Xy) dy. 

Note that the coefficient fi does not have any effect on a linear operator. 

We will call L\ to the largest scale invariant class contained in the class £* defined in (|2.6| - (|2.7D . 
This is the class of integro-differcntial operators with kernels K such that 

^--)^<K(y)<{2~a)^ (5.1) 

\VK(y)\ < C\y\~ n ~' T - 1 in R" \ {0}. (5.2) 

From Theorem 12. 7\ an equation Iq(u,x) = has interior C 1,a estimates for some a > if Io is 
uniformly elliptic respect to the class C\. 

Our main theorem in this section states essentially that if an equation 1^ (u, x) = is uniformly 
elliptic respect to a scale invariant class with interior C 1 '" estimates and we have another equation 
I(it, x) = f(x) for which the operators I(— , x) stay close to T ', then this equation also has interior 
C 1,a estimates for any a < min(a, a — 1). Since our proof involves successive rescaling, we need to 
measure closeness at every scale, so we define a distance in between operators that takes scaling 
of order a into account. 

Definition 5.1. Given a £ (0,2) and an operator I, we define the rescaled operator as above 

I„>,i) - AV^'MA" 1 -), Xx) 
The norm of scale a is defined as 

I« _t(2) 



sup 

A<1 



t(1) _ t(2) 

h.x H.x 

:i 

where \\.\\ is the norm defined in Definition \2.8l 

The purpose of the rescaled operator is to stress that if u solves the equation l(u,x) = f{x) 
in B\, then the rescaled function w(x) = fiu(Xx) solves an equation of the same ellipticity type 
l^x(w,x)=X°f(Xx) inBi. 

The following is the main theorem of this work. 

Theorem 5.2. Assume a > o~o > 1. Let be a fixed translation invariant nonlocal operator in 
a class C C £o(o~) with scale a and interior C 1 '" estimates (for example C — C\). 
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Let I^ 1 ) and I^ 2 * 1 be two nonlocal operators, elliptic with respect to Cq{<j), and assume that 



l(0) _ j(i) 



< n 



for some rj > small enough and j = 1,2. 

Let u be a bounded function that solves the equation 

l {l) {u,x) >fi(x) inB x 
I {2) (u, x) < f 2 (x) inBx 



for a couple of bounded functions f\ and f 2 . 

Then u £ C 1,a (Bi/2) for any a < min(a, ctq — 1) and 



l«llcn.-(fl va ) < C (\\ u \\l° 



> ) + \\M\l~ { b,) + \\M\l~(bJ 
and the estimate depends only on a , X, A and dimension, but not on a. 

The proof of this theorem is based on an approximation argument inspired in [3J using Lemma 

ETT1 

Proof. By scaling the problem, we can assume without loss of generality that 



ll/llli. 




< »?, 


II/2IL0 




< »/, 


IMLo 




< 1, 



and m solves the equation in some large ball B^r. 

We will show that there is a A > and a sequence of linear functions 



such that 



l k (x) = a k + b k ■ x, 



sup|u-Z fc | < X k ^ 1+a) 
s Afc 

|ofc+i-o fc | < X k(1+a) 
X k \b k+1 -b k \<C 2 X k ^ 



(5.3) 

(5.4) 
(5.5) 



It is a standard computation that these conditions imply that u is C 1 '" at the origin (see [3]). 
Let Iq = 0. We proceed by induction. Assume we have (|5.3[) - (|5.4[) - (|5.5[) up to some value of k. 
We will now show it for k + 1 . Consider 



w k (x) = 



[u-l k ](X k x) 



A fc(l+a) ' 

then since C has scale a, w k solves an equation of the same ellipticity type 

I^K.x) = I^} H1+ahxk (w k ,x) > X k ^- x -^h{X k x) (and < for if and f 2 ) 
Since a — 1 — a > <7q — 1 — a > 0, the right hand side becomes smaller as k increases. Moreover, 



I 



U) _ T (o) 



k 



< 



lU) _ 1(0) 



< V- 
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By the inductive hypothesis, we have that \wk\ < 1 in B\. Let a < ol\ < min(a:,<7o)7 we will 
show that we can construct the sequence Ik so that we also have 

w k (x) < \x\ 1+ai if x$Bi 

We choose R from Lemma 14.71 so that we have that Wk is Holder continuous in Br \ B\ and 
we can apply Lemma 14.71 to the function h that solves 

l[ 0) h = in B x 

h(x)=w k (x) inR n \B 1 

then \wk — h\ < e(rj,R) in B\. 

Since Ij, is translation invariant and elliptic with respect to C, it has interior C 1 ' estimates. 
Let I = a + b-x be the linear part of h at the origin. Since \h\ < l+e(r],R) in Bx, \a\ < l+e(r], R). 
By the C 1,a estimates, b < C 2 . We also have that \h(x) — l(x)\ < C3|x| 1+Q in Bx/2- So, we have 
the following estimates: 

\w k (x) - l(x)\ < e( V ,R) + C 3 A 1+a in B 1/2 (5.6) 

\w k (x)-T(x)\ < e(?7, i?) + 2 + C 2 mBi\B 1/3 (5.7) 

\w k (x)-T(x)\ < e(r], R) + 1 + \x\ 1+ai + C 2 \x\ in K" \ B\ (5.8) 

We will choose A small to be determined below, and then rj and R so that e(i], R) < A 1+Q . The 
question is how to pick A so that we can complete the inductive step. We have 

i k+1 (x) = i k (x)+x k ^ +a n^), 

(u-l k+ i)(X k x) _ (w k -l)(\x) 
Wk+1 - A fc ( 1+C ") ~ A!+ a 

Let us analyze the bounds (|5.6[) . (|5.7|) and (|5.8p in the rescaled setting. We get 
\w k +i(x)\ < A a - Q + C 3 A a - Q in B x - 1/2 

\w k+ i(x)\ < \ & - a + (2 + C 2 )\~ {1+a) in Bx \ B x -i /2 

\w k +i(x)\ < \ & - a + 1 + |A| Q1 - Q |x| 1+Ql + C 2 |Ar Q |x| in W l \ B A _i 

We start by choosing A small enough so that (1 + C?,)\ a ~ a < 1 and we have \wk+i | < 1 in Bx- 
If we choose A small enough, we also have the following estimates in each of the three rings 

\wk+i(x)\ < \ & - a + C 3 X & - a < 1 < \x\ 1+ai in B x - 1/2 \ Bx 

\w k +i{x)\ < A 3 "" + (2 + C 2 )A-( 1+Q ) < c\ ai - a \- (1+ai '> < \x\ 1+ai in B x -i\B x -i /2 
\w k+ x(x)\ < A a - a + l + |A| ai - a | a; | 1+ai +C 2 |A|- Q |x| < |x| 1+Ql ■mM n \B x -i 

so that |wfe+i| < outside Bx- Moreover we observe that 

\a k+ i - a k \ < X k(1+a) 
\ k \bk+i-bk\<C 2 \ k( ~ 1+ ^ 

where C 2 is the same as above and depends only on the C 1,a estimates. We also have u(x) — 
l k+ x(x) = A fe ( 1+Q > - thus 

\u(x) - l k+ x(x)\ < A( fe+1 )( 1+a ) in B xk+1 . 
This completes the inductive step and the proof. □ 
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6 Applications 



Since the main theorem in scction[5]was stated in its more general way which is somewhat abstract, 
we provide concrete applications in this section aimed at showing the strength of this result. 

6.1 Linear equations with variable coefficients 

We start by stating the result of the theorem in the simplest case when the equation is linear and 
close to an operator in C\. This is the integral version of the classical Cordes-Nirenberg theorem. 

Theorem 6.1. Let a > 1 and let u be a bounded function that solves the equation 

x) := J 5u(x, y)(2 - dy = f(x) m B t 

E" 

Assume that f G L°°(]R n ) and \a(x,y) — ao(y)\ < V for every x G B\ and some small r\ > 0. 
Assume do is a bounded function such that K(y) = (2 — a)a^{y) /\y\ n+a satisfies (|5.1[) and (|5.2[) . 
Then for any a < a — 1, u € C 1,a (B 1 / 2 ) if f] is small enough and 

IMI C i,«(b 1/2 ) < c(.\W\\l°°(r») + II/IIl^b!))- 

Proof. We apply Theorem 15.21 In this case is given by 

i(°) (u ,,) = /M„)^Pd, 

By assumption, this operator I^ -* belongs to C\, which is a scale invariant class with interior 
C 1,a estimates. In this case it is easy to see that since the equation is linear and the coefficients 
do not depend on x, the derivatives of the solution solve the same equation and the solutions are 
actually C 2 ' a . So in particular has interior C 1 ' 1 estimates. 

The weight u> must be a weight that controls the tails of the kernel. In this case a natural 
choice is 

1 

w(y) 



(i + M) 



n+cr 



Now we estimate III — I 



!_!(0) 


= sup 







Il,A - (I (0) )l,A 



We compute 



(Ii,x-(l (0) )i,x)(u,x)= 5u(x,y) 



(2 - a)(a(Xx,Xy) - ao(Xy)) 

\y\ n+a 



with |a(Aa;,Ay) — 1| < ry. From Definition 

Il,A-(I (0) )l,A 



= sup 

u as in Definition l2"2l 

1 

— SU P 77 

u as in Definition [O] 1 + M 



l hX {u,x)-(l^) hX (u,x) 



1 + M 

\Hx,v)\ I^IJ 1 dy 



(6.1) 
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In Definition we take functions it such that ||u|| £ i/ R „ w s and \u(y) — u(x) — (y — x) ■ Wu(x)\ < 
M\x — y\ 2 for y G Note that these two conditions together imply \u(x)\ < CM for some 

universal constant C. 

From \u{y) — u{x) — (y — x) ■ Vit(x)| < M\x — y\ 2 , we obtain 

< CM 



On the other hand, from ||it|| L i( Rn u \ and \u(x)\ < CM we bound the tail of the integral 
\5u(x,y)\^±dy< J \u(x + y)\ + \u(x - y)\ + \u(x)\^0 dy 

»"\Bi R"\Bi 

< CM 

Replacing in (|6.1[) . Ii.a — Ii°a < C ? ? f° r anv thus ||l — < Cr\. So if we choose r\ 
small enough we can apply Theorem 15.21 and conclude that the equation I(u,x) = f has interior 
C 1 '" estimates for any a < a — 1. □ 

It is interesting to note in the theorem above that if we fix an a > 0, the constant C in Theorem 
16. II does not blow up as a — > oo. 

6.2 Nonlinear equations with variable coefficients 

By combining Theorem 15 . 2 1 with Theorem 12 . 71 we obtain the following result 



Theorem 6.2. Let a > <jq > 1 and let £i(<r) be the class defined in (|5.ip - ()5.2|) . 

Assume I*- -* is translation invariant nonlocal operator uniformly elliptic respect to C\{a) and 
I uniformly elliptic respect to Cq. Assume ||l(— ,x) — I*- -^^ < i) for every x £ B\ and some small 
T) > 0. Let f G L°°(Bi) and u be a bounded function that solves the equation 

I(it, x) = /(x) m £?i. 

TTien it G C 1,a (B 1 / 2 ) for some small a > and 

Mci.«(b 1/2 ) < C(\\u\\ LOO(Mn) + ||/|| L= o (Sl) ). 

and the estimate depends only on o~o, A, A, C\ and dimension, but not on a. 

A concrete example of a nonlinear operator I such that Theorem 16.21 applies would be 



I(n,x) = inf sup / Mx + y) + u(x - y) - 2u{x)) ^ ~ -)My) +aM^ dy 

R" 

such that we have for some small rj, 

\a a p(x,y)\ < T] for every a,0 
A < ao(y) < A 
|V ao (2/)| < CIj/I" 1 

Of course, a regularity result that requires coefficients not to vary much with respect to x 
implies a regularity estimate every time the coefficients are continuous in x. We state it in the 
following theorem. 
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Theorem 6.3. Given a > <7q > 1, let l(u, x) be given by 

l(u,x) =infsup / 5u(x,y) — — %^r^ ^ dy 
a p j \y\ n+rT 

such that for every a, j3 we have A < a a p(x,y) < A and V v a a p(x,y) < Ci/((2 — er)|j/|) and also 
\aap(xi,y) — a a p{x2,y)\ < c(\x± — X2I) for some uniform modulus of continuity c. Then solutions 
u to the equation 

l(u, x) = f in 
are of class C 1,a (B 1 / 2 ) for some small a > and 

Mcl."(B 1/2 ) < C(l|w|ll,-(Kn) + ||/||i=o (Bl) ). 

and the estimate depends only on o~o, X, A, C\, the modulus of continuity c and dimension, but 
not on a . 

Proof. Around each point xq S -B1/2 we can find a ball B r (xo) so that |a Q| g(x, y) — a a /3(xo,y)\ < 
c(r) < rj. This implies that ||I(— ,x) — I(—,Xo)\\- < Ctj as in the proof of Theorem 16.11 and we 
apply Theorem 16.21 with I^ ' = I(— , xo), scaled in B r (x). □ 

6.3 Nonlinear equations with constant coefficients but non-differentiable 
kernels 



It is interesting to note that Theorem 15.21 is useful to obtain new results even in the translation 
invariant case. In Theorem 12. 7\ it is required that every kernel must be differentiable away from 
the origin. That assumption can be relaxed in the following way. We can actually obtain a C 1 '" 
estimate for nonlocal equations that are uniformly elliptic respect to the class C given by operators 
with kernel K such that 

A < ai(y) < A 

Nl < v 

|Voi(y)|<^ inM»\{0} 

Theorem 6.4. Assume a > <7o > 1. Iff] is small enough (depending on A, A, C\ and dimension) 
then the nonlocal equation 



l(u,x) = mf sup / (u(x + y) + u(x - y) - 2u(x))K a p(y) dy = f(x) 
a 13 J 

R™ 

has interior C ,a estimates if all kernels belong to the class C above. 

Proof. Let L £ C be an operator with kernel K, we can write it as K = K 1 + K 2 where A' 1 = 
(2 — a) |°[n+l an d correspondingly L = L 1 + L? . But then we see that ||i — L x \ a < crj- So we 
write 

I^ (u, x) = inf sup L^gu(x) 

a p 

and we have ||I — || ct < ci] so that we can apply Theorem 15.21 □ 

This theorem applies for a class that is still smaller than Cq. It would be interesting to 
determine whether the class £0 has interior C 1 '" estimates or not. We leave that question open. 



23 



6.4 Nonlinear equations close to the fractional Laplacian. Very regular 
solutions. 



We can obtain another new result in the translation invariant case from Thcorem l5.2[ In this case 
we look at non linear equations that are sufficiently close to the fractional Laplacian and we obtain 
C 2 ' a regularity. We provide a theorem on the regularity of fully nonlinear translation invariant 
nonlocal equations when the ellipticity constants arc sufficiently close to each other. This is an 
improvement on Theorem 12.71 when these conditions arc satisfied. 

Theorem 6.5. Assume a > o~o > 1. There is a r\ > and po > so that if A < I + r\, A > I — r\, 
I is a nonlocal translation invariant uniformly elliptic operator with respect to C* (given in (|2.6p - 
(|2.7p ) and u is a continuous function in B\ such that u G L 1 (M n ,uj), lu = in B\, then there is 
a universal (depends only on n and oq) a > such that u G C 2Jra (B\j2) and 

u c 2 +"(B 1/2 ) < C (supM + \\u\\ Ll(Rn uj) + |/0 

for some constant C > (where by 10 we mean the value we obtain when we apply I to the function 
that is constant equal to zero). The constant C depends on ao, n and the constant in (|2.7[) . 

Proof. From Theorem 12.71 we know that u G C 1 ' a (B 3 / 4: ). In particular the function is differ- 
entiable. Let w ~ u e be a directional derivative. We proceed as in the proof of Theorem 12.71 
(Theorem 12.1 in [2]) and write w — w\ + W2 where wi vanishes in -B5/8 and w\ solves 

M^wi > -C in £5/8 

M^wi < C in £5/8 

The difference with the proof of Theorem 12.1 in [2] is that now we apply Theorem 15.21 in- 
stead of Theorem 12.61 Since A > 1 — T) and A < 1 + 77, then ||M£ t<>i + (-A) CT || < crj and 
||M^ wi + (— A) CT || < or). So Theorem 15.21 tells us that w = u e is C 1,a in S1/2, and thus 
ueC 2 > a {B 1/2 ). □ 



Theorem 15.21 is also useful to reduce the equation to its leading order part. A simple example 
would be the following. 

Theorem 6.6. Let a G (1, 2) and r > 0. Let u be a bounded function that solves the equation 
l(u,x) ^ -(-A)"/ 2 u{x) + mi sup I fa(x,y)(2-cr) |y| , r ° a f (y) dy = in B x 

a p J \y\ n+,T 

R™ 

Assume that f G L°°(R Tl ), r > 0, \a(y)\ < C, and \Va(y)\ < C/\y\. Then u G C 2 ' a (B 1/2 ) for 
some a > and 

ll U llc 2 ' = (Bi /2 ) - C II M IIl°°(R") ■ 

Proof. First of all let us notice that we can place the fractional Laplacian inside the inf sup and 
rewrite the equation as 

1^) = inf sup f 8u{x, V ){2-a) C ° + ] y ll a / {y) dy = Q in B 1 
a p J \y\ n+ ' J 

where c a is the normalization constant of the fractional Laplacian which remains bounded below 
and above for a G (1,2). 
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When wc make a rescaling of the equation of order a we obtain 

Ii, A =mfsup / bu (x, y) 2 - a) j-r— — dy. 

Then, for A small enough, we will have ||li.A + ( — ^0 || < V an< i we proceed as in the proof of 
Theorem 16.51 □ 

Appendix 

In this appendix, we prove Lemma 13. II We can give an intuitive proof avoiding lengthy computa- 
tions by using Lemma 14.31 We notice that we did not use Lemma |3. II to prove Lemma 14.31 

Proof of Lemma \3.1\ For some r > 0, we will show that if we consider u a (x) = ((|x| — 1) + ) Q and 
xq = (1 + r)ei, then M + u(xq) < if a is small enough. By rotation this implies the result for 
any x such that \x\ = 1 + r. Ifl<|x|<l + r, we can still get the result by a scaling centered at 
e\ (i.e. by considering u(x) = {{\x + Aei| — (1 + A))+) Q for A such that \x \ — 1 = r(l + A)), since 
that scaling would only decrease the values of the function elsewhere. This observation implies 
also that if M + u(xo) > 0, then also M + u(x) > for all x such that \x\ < \xq\. 
Let M be a large number and let Vj be 

, , f u a (x) — 1 \ 

v a (x) = max , — M . 

V a ) 

The key observation is the following 

lim v a {x) = max(log(|x| — 1), — M) = Im{x) 

ol — ^0 

and this convergence holds uniformly on compact sets. 

As a —> 2, M+ converges to a uniformly elliptic second order differential operator that is 
comparable to the maximal Pucci operator with ellipticity constants depending on A, A and 
dimension, that we call ■ It is a simple computation that M + (D 2 log(|x| — 1)) —> — oo as 
\x\ —> 1 (where M + (D 2 u) stands for the standard Pucci operator). So if we pick r small enough 
and M larger than — logr, then there will be a g\ such that M+Zm < — 1 if <r > oi. Moreover, 
by choosing M large enough, we can also make sure that M^Im < — 1 for a < o\. This is how we 
pick r and M, 

Assume that the result of the lemma was not true. So there would be a sequence ctj — ¥ and 
Gj G [ct ,2] such that M+.m Qj . (x ) > 0, and thus M+v aj > in the annulus B 2 \ Xo \ \ B\ Xo \. We 
extract a subsequence o~j k that converges to some number a € [ctq > 2] . We apply Lemma 14.31 and 
we obtain that M^Im(xq) > 0. But this is a contradiction since M£Im{xq) < — 1. □ 
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